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Abstract. We derive practical expressions for a rapid and 
accurate evaluation of electric and thermal conductivities 
and thermopower of degenerate relativistic electrons along 
quantizing magnetic fields in outer neutron star crusts. We 
consider the electron Coulomb scattering on ions in liquid 
matter, as well as on high-temperature phonons or charged 
impurities in solid matter. We propose also a reasonable 
semi-quantitative treatment of low-temperature phonons. 
The transport properties are expressed through the energy 
dependent effective electron relaxation time averaged over 
electron energies. We calculate this relaxation time, using 
the theoretical formalism of the previous work, obtain ac- 
curate fitting expressions, and propose an efficient energy 
averaging procedure. We create a computer code which 
calculates the longitudinal transport properties of degen- 
erate electrons in strong magnetic fields for any parame- 
ters of dense stellar matter of practical interest. We anal- 
yse quantum oscillations of the transport coefficients ver- 
sus density at various temperatures and magnetic fields. 
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1. Introduction 

Transport properties of neutron star crusts are important 
for studies of thermal evolution (cooling) of neutron stars, 
evolution of their magnetic fields, etc. For instance, the 
heat in the outer crusts of magnetized neutron stars is 
mainly transported by electrons along the magnetic fields. 
As a rule, the electrons in the crust are strongly degener- 
ate, and they may be relativistic. The main electron scat- 
tering mechanisms are the Coulomb scattering on ions in 
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the liquid phase, and the scattering on phonons or charged 
impurities in the solid phase; the magnetic field can quan- 
tize electron motion. 

In the present work we study the electron transport 
properties in the outer crusts of neutron stars along quan- 
tizing magnetic fields. The problem has been considered 
earlier by several authors. The first were Canuto and 
his colleagues in 1970s; the results were summarized by 
Canuto & Ventura (1977). An adequate kinetic equation 
for the electron distribution function was proposed by 
Yakovlev (1980), and used for calculating the transport 
coefficients by Yakovlev (1980, 1984), Hernquist (1984), 
Van Riper (1988) and Schaaf (1988). Although much work 
has been done, these studies are not fully complete. First, 
the approach based on the distribution function formalism 
is not invariant: the results depend slightly on the choice 
of the electron wave functions (due to the spin degener- 
acy) in the magnetic field. Second, the results are not easy 
for practical use. Third, the previous studies neglected 
the Debye- Waller reduction of electron-phonon scattering 
rate in solid magnetized matter, although the importance 
of the Debye-Waller factor was demonstrated by Itoh et 
al. (1984, 1993) for the non-magnetic case. 

Recently Potekhin (1995, hereafter Paper I) has devel- 
oped an invariant formalism of the longitudinal electron 
conduction problem based on the spin polarization den- 
sity matrix of the electrons. The results of Paper I have 
been compared with those obtained using the traditional 
methods with various electron basic functions. The lat- 
ter methods are simpler for calculations, while the den- 
sity matrix formalism enables one to choose which basis 
is more adequate at given parameters of stellar matter. 
Paper I has presented also efficient numerical techniques 
for calculating the electron transport coefficients includ- 
ing the Debye-Waller factor and demonstrated that this 
factor is much more significant when the field is strongly 
quantizing. 
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In this paper we obtain practical formulae for a simple 

and rapid evaluation of the longitudinal transport coeffi- 
cients using the theory developed in Paper I. 

2. Physical conditions 

We consider degenerate layers of an outer neutron star 
crust (at densities p ^ 4 x 10^^ g cm"'', below the neu- 
tron drip). Matter of these layers consists of electrons and 
ions. We study not too low densities (see below), at which 
the electrons are nearly free, and the ionization is com- 
plete due to the high electron pressure. The ions give the 
major contribution into the density while the electrons - 
into the pressure. For simplicity, we mainly consider one 
component plasma of ions (=nuclei {A, Z)); the electron 
and ion number densities are related as rie = Zm. A com- 
prehensive study of thermodynamic properties of magne- 
tized neutron star crusts has been performed recently by 
Rognvaldsson et al. (1993) using the Thomas-Fermi ap- 
proximation. Since we avoid low densities, we base our 
consideration on the approximation of free electrons. This 
will enable us to evaluate the transport properties of mat- 
ter. We will include the effects of non-ideality of electrons 
in a phenomenological manner in Sects. 5 and 6. 

Let n be the electron chemical potential (including the 
electron rest-mass energy, m,c?). The domain of strongly 
degenerate electrons we are interested in corresponds to 
T <^ Tp, where Tp = (/U — mc^)/fcB is the degener- 
acy temperature and fee is the Boltzmann constant. If 
{fi — mc^) <C mc^, then the electron gas is non-relativistic, 
while for fj, » mc^ it becomes ultra-relativistic. 

In the absence of the magnetic field, one has the famil- 
iar result 




B = 

B=lO\l G 
13 B=10^^ G - 
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Fig. 1. p - T diagram of '""^Fc matter for several magnetic fields 
strengths. Tp - electron degeneracy temperature, Ti - ion cou- 
pling temperature (F = 1), Tm - melting temperature, Tp - 
ion plasma temperature, Tb ~ electron quantization tempera- 
ture; arrows indicate strongly quantizing densities ps. Lines A'^ 
restrict the low-T low-p domain of incomplete ionization and 
electron gas non- ideality. Tp, Tb, and N curves depend on the 
magnetic field. Dashes show the curves for B = 0, solid lines 
for B = 10^^ G, and thick long-dash lines for B = 10^^ G. 
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is the Fermi-Dirac distribution function, and e = 
c^/ (mc)2 +p'^ is the electron energy. In the case of strong 
degeneracy, one has /U « /xq, where 



Mo 



= ^m^c"^ -h c^pIq, Pfo = ?i(37r^ne) 



1/3 



(3) 



and Pfo is the field-free electron Fermi momentum 
(pfo/(wc) « 1.009(p6.^/^)^/^, P6 being density in units 
of 10^ g cm"''). The appropriate degeneracy temperature 
is plotted in Fig. 1 by short-dashed line. For small p and 
T, the effects of incomplete ionization and electron gas 
non-ideality become important. This domain is shown 
schematically by the short-dashed line N. The line cor- 
responds to s = \sa\, where e is the mean energy per 



electron in the free electron gas, and £„ is the mean en- 
ergy per electron for isolated atoms in the Thomas-Fermi 
approximation (e.g.. Landau & Lifshitz 1976). 

Let the magnetic field B be directed along the ^;-axis. 
We shall use the Landau gauge of the vector potential: 
A = {—By, 0, 0). Quantum states of a free electron can be 
labelled by four quantum numbers {px,Pz,n, s), where p^ 
determines the t/-coordinate of the electron guiding center, 
Vb = Px/imuJB), Pz is the longitudinal electron momen- 
tum, n=0,l, 2,. . .enumerates the Landau levels, and s is 
a spin variable. Here cob = eB/{mc) is the electron cy- 
clotron frequency, and (— e) is the electron charge. The 
ground Landau level n = is non-degenerate with re- 
spect to spin while the excited levels n > are doubly 
degenerate. Various sets of electron basic wave functions 
are analysed, e.g., in Paper 1. An energy e of a relativistic 
electron in the magnetic field is 



Emc^ 



where E is the natural dimensionless energy. 
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The number density of free electrons in a magnetic 
field is 



He = 



mtOB 
(2^ 



/+00 
-oo ^ „ 



(5) 
(6) 



s n— 



Here Pnis) = mcPn{E) is the largest z- momentum of an 
electron with energy £ at a given Landau level, n{E) = 
Int(i^) is the highest Landau level occupied by an electron 
with energy s, 



Pr,{E) = VE^-l-2bn, u = {E^-l)/{2b), 



(7) 



Int(j^) is integer part of v, and b is the dimcnsionless 
magnetic field expressed in units of the 'critical field' Be 
{huBc = mc^), 



B 



2 3 

6 = ^, Be = = 4.AU X 10^3 G. 

Br en 



(8) 



Mass density p equals p w mini w {A/Z)mune, where 
mi is the ion mass, rriu = 1.660 x 10~^^ g is the atomic 
mass unit, and A/Z is the number of baryons per electron. 
Equation (5) determines TIq clS 3/ function of T, jjL and B. 

Simple analysis of Eq. (5) yields the following. Strongly 
degenerate electrons populate the only ground Landau 
level when density is not too high, p < pB, where pb = 
2.07 X 10^ b^^'^{A/Z) g cm-3. In this case the Fermi mo- 
mentum is 



Pf = 



^ fPFO\- 

= tjtPfo 



(9) 



In the domain of p <C pb the magnetic field strongly 
reduces the electron Fermi energy and widens the re- 
gion of incomplete ionization and electron gas non-ideality 
(Yakovlev 1984, Hernquist 1984, Van Riper 1988), as 
shown in Fig. 1 by solid and long-dashed lines. In this 
case, lines N correspond to e = \£a\, where depends on 
B (Rognvaldsson et al. 1993). 

If B is fixed, then p, increases with p, and degenerate 
electrons populate higher Landau levels n. If p ^ ps, 
many Landau levels are populated, and the Fermi energy 
is almost independent of B (pp ~ Pfo)- 

It is convenient to introduce the temperature 



Tr = 



kB 



1.34 X 10^Si2 



mc 



K, 



(10) 



where ujg = eBc/e* is the gyrofrequency of an elec- 
tron with the typical energy e* = max(£, /x), and B12 = 
B/(10^^ G). When T :$> Tb the electrons occupy many 
Landau levels for any p due to high thermal energy. In 



this case the thermal widths of the Landau energy lev- 
els (^ fceT) are higher than the inter-level spacing, and 
the magnetic field acts as non- quantizing, regardless of the 
electron degeneracy. 

In the domain of T ^ Tb and p < pb (separated into 
the subdomains of degenerate and non-degenerate electron 
gas) the electrons populate mostly the ground Landau 
level. In this domain, the magnetic field is strongly quan- 
tizing, and it affects significantly all properties of matter 
(e.g.. Van Riper 1988, Yakovlev & Kaminker 1994). 

Finally, in the domain of T ^ Tb and p^ pb the elec- 
trons are degenerate and populate many Landau levels but 
the inter-level spacing exceeds k^T. Then the magnetic 
field is weakly quantizing. It does not affect noticeably the 
bulk properties (pressure, electron chemical potential) de- 
termined by all the electron Fermi sea but it affects the 
transport properties determined by thermal electrons near 
the Fermi level. 

The state of ions can be characterized by the ion cou- 
pling parameter 



r = 



ak^T 



^2 /p6\i/3 



(11) 



where a = [3/(47rni)]^/^ is the ion sphere radius, and T7 
is temperature in units of 10^ K. 

At high T, when F <C 1, the ions form a classical Boltz- 
mann gas. With decreasing T, the gas gradually (without 
any phase transition) becomes a Coulomb liquid. The liq- 
uid is formed (Hansen 1973) at F « 1, i.e., at T k, T\ 
(Fig. 1). The liquid solidifies into the Coulomb crystal 
(Nagara et al. 1987) at F = 172 {T = T^). According 
to Fig. 1, we always have strongly coupled Fe ions, if elec- 
trons are degenerate {T\ > Tp). However, for light ions, T\ 
can be lower than Tp. 

At low T, zero-point ion vibrations become important 
in ion motion. These effects are especially pronounced if 



T < Tp, where 



Tp = Uw^/kB « 7.832 x \{f{Z/A)pl''^ K 



(12) 



and LOp = ^ AirZ'^e^ni/mi is the ion plasma frequency. 
The amplitude of zero-point vibrations is commonly much 
smaller than the typical inter-ion distance, a. With in- 
creasing p, the amplitude-to-a ratio becomes larger, so 
that the vibrations can prevent crystallization at high p 
(Mochkovitch & Hansen 1979, Ceperley & Alder 1980). 
This effect is especially pronounced for H and He. Note 
that the Debye temperature of the Coulomb bcc crystal is 
Td = 0.45Tp (Carr 1961). 

We assume that the magnetic fields do not affect the 
properties of the ion plasma component, for instance, the 
phonon spectrum of the crystal. This is so (e.g., Usov et al. 
1980) if lob = ZeB/{mic) < Wp, i.e., if B «: lO"^ G. 

3. Transport properties 
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3.1. Electron scattering mechanisms 

The electron transport properties are determined by the 
electron scattering mechanisms. Consider three important 
cases when the electron scattering is almost elastic (an 
energy transfer in a colHsion event is <C k-QT): (i) the 
Coulomb scattering on ions in the liquid or gaseous phase 
(T > Tni, Fig. 1), (ii) the scattering on high-temperature 
phonons (Td ^ T < Tj„); (iii) the Coulomb scattering on 
charged impurities in the lattice (important for T <^ Td). 
The impurities represent ions of charge Zi^p ^ Z im- 
mersed accidentally in lattice sites. Accordingly, our re- 
sults cover a wide range of temperatures. 

The electron-ion scattering can approximately be de- 
scribed (e.g., Yakovlev 1984) by the Debye-screened 
Coulomb potential. Its Fourier image U{q) is given by 



\Uion{q) 



(9: 



ion \ 2 



2 I 2 



(13) 



where qs is an effective screening wavenumber (inverse 
screening length), q{ and qc are, respectively, the ion and 
electron screening wavenumbers. According to Yakovlev 
& Urpin (1980) and Yakovlev (1984), 



9i 



(2a/3)^ 



' D' 



where is the Debye ion radius, ~ Vi/i^p, Vi = 
^ k-QT/rrii is the thermal ion velocity. The Debye radius 
dominates [qi ~ I/'^d) in the gaseous regime T » Ti, 
while qi w 1.5/ a in the most important case of strongly 
coupled ion liquid (Tm < T <, Ti). Equation (|lj) allows 
us to reproduce the B = transport properties using the 
model potential ( p^ instead of the exact plasma-screened 
Coulomb potential. 

For the scattering on high-temperature phonons, one 
has (e.g.. Paper I) 



2 _ r ^wZe^ 
V 1 



exp [-2W{q)] 



(15) 



where = u_2a^/r is the mean squared thermal dis- 
placement of ions from their lattice sites, U-2 is a numer- 
ical factor determined by the phonon spectrum (u_2 = 
13.0 for the body-centered cubic (bcc) lattice, see, e.g., 
Mochkovitch & Hansen 1979), e"^^ is the Debye- Waller 
factor, 2W{q) « (rTg)^/3. Theoretical formalism for in- 
corporating this factor into the transport properties of 
magnetized matter was developed in Paper I. Note that 
Eq. (^ is vahd for Td T < T^. 

The scattering on phonons remains to be the dom- 
inant electron scattering mechanism at lower tempera- 
tures, T <C Tq. The low-T phonon scattering is inelastic 
(electron energy transfer in a collision is ^ k^T), and, 
strictly speaking, it cannot be described using the formal- 
ism of Paper I. However this scattering can approximately 
(semi-quantitatively) be treated as elastic (Yakovlev & 



Urpin 1980, Raikh & Yakovlev 1982) as long as the Umk- 
lapp processes dominate over the normal processes of 
electron-phonon interaction. The Umklapp processes are 
known to be most important at high temperatures, un- 
der typical conditions in neutron star crusts. However 
they are "frozen out" , and the electron scattering be- 
comes essentially inelastic. For instance, for nonquantiz- 
ing fields, this happens at T < Tf/ ~ TpZ'^/^e'^ / (hvp) « 
[Z2/3 + 0.98{A/ pe)^/Y/'^(Tp/l37), where vf is the elec- 
tron Fermi velocity. We propose to extend the formal- 
ism of high-temperature phonons to T <C Td and obtain, 
thus, the reliable estimates of low-temperature transport 
properties at Ty ^ T ^ Td. For this purpose, as can be 
shown, for instance, from the results of Raikh & Yakovlev 
(1982) and Baiko & Yakovlev (1995), the scattering po- 
tential ( |l5| ) should be modified in two ways. First, the pre- 



exponent factor r^ should be replaced by 



Tl 



■G(t), 



where G{t) = t/ i/P + i^ describes the reduction of ther- 
mal ion displacements at low T, t = T/Tp, and to ^ 0.1 
is a numerical parameter {to ~ 0.132 for the bcc lattice). 
Second, one should remind that the Debye- Waller factor 
is determined by the total (thermal -I- zero point) mean 
squared ion displacement which reduces to the ther- 
mal squared displacement at high T. Thus we should 



(14) replace by r. 



2 

T2 



^,{1 + exp{-9.lt)[u-i/{2u2t)]} in 



the Debye-Waller exponent (Baiko & Yakovlev 1995) for 
it to be accurate at low T. In this case, u-i is another 
parameter of the phonon spectrum (u_i = 2.800, for the 
bcc lattice). The above modifications do not violate the 
formalism of Paper I. 

Finally, the scattering on impurities corresponds to 
(e.g., Yakovlev & Urpin 1980) 



|C/imp(q) 



4^(Zi„,p - Z)e^ 



(16) 



In this case the screening wavenumber is {qf^^Y = q^ + 
(7?„p , where qc is given by Eq. ( p^ ) , and qimp is an inverse 
impurity correlation length. The scattering on impurities 
is similar to that on ions. This scattering acts at T < Tm, 
just as the scattering on phonons, but actually it domi- 
nates at very low T (see Yakovlev & Urpin 1980). 



3.2. Transport coefficients 

Let dfji/dz and dT/dz be, respectively, weak and locally 
constant electric field, electron chemical potential gradi- 
ent, and temperature gradient along B. They induce the 
electron electric and thermal currents with the current 
densities J and Q: 



J 



(7[8 + - — 

e oz 



dT 

oz 



~f3T 



e oz 



A— 

dz 



(17) 



Here a is the longitudinal electric conductivity, while /3 
and A are two other auxiliary transport coefhcients. The 
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appearance of f3 in the expressions for J' and Q reflects 
the Onsager symmetry principle. For practical use, the 
expressions (17) can be rewritten as 



^ Idfi J dT 
e oz a oz 



Q 



dT 

-aTJ - K—, 
oz 



where 



P/a, X-T{l3^/a) 



(18) 



(19) 



are, respectively, the longitudinal thermopower, and ther- 
mal conductivity. The transport coefficients a, (3, and A 
are convenient for calculation (see Eq. (22) below) while 
a, a and n directly enter the equations which govern the 
distributions of temperature and magnetic field in neutron 
stars (e.g., Urpin & Yakovlev 1980). The coefficients cr, k, 
and a determine fully the electron transport of charge and 
heat along B. 

Let Pns'si^iPz) be the spin polarization density matrix 
(Paper I) of electrons for the case when the electron gas is 
slightly non- uniform along B. Here s and s' are the spin 
variables (Sect. 2). In the linear regime, deviations from 
the equilibrium are small. The zero-order density matrix 
is f{e)Sss', where /(e) is a local Fermi-Dirac distribution 
(H) which depends on z parametrically through /i and T. 
In the ffist approximation, according to Paper I, 



Pnss'{z,Pz) = f{e)5ss'+r]^ 
e~ PL dT 



dfie) 



T dz 



de 



c dp 

OZ 



(20) 



where rj =sig\i{pz), and the functions iprinss'i^) determine 
non-equilibrium corrections to the density matrix; / is an 
electron scattering length: 



mc^ huj B 



ph 



27miZ2e4 ' 
27mimp(^imp - ZYe^' 



he 



(21) 



''imp 



being the number density of impurities. The set of 
equations for iprjnss'i^) has been derived in Paper I. Pa- 
per I presents also the mathematical formalism for solving 
these equations. After introducing the scale lengths (21), 
the equations describing the scattering on ions and impu- 
rities appear to be formally the same, i.e., there is no need 
to consider these scatterings separately. We shall refer to 
the scattering on ions and impurities as the Coulomb (C) 
scattering. If the functions ifjjnss' (^) are found, the trans- 
port coefficients in Eqs. (17) can be calculated as (Paper I) 



'cr 



'™^/ \eie~p}/T 

l^'ncA ie- pf/T 



(27rn)2 



X *(£)(-|7 )de. 



(22) 



where 



n(E) 



= ^ E E E ^''-(^)- (23) 

r)=±l n=0 s=±l 

It is also convenient to introduce the function ^'(£'), 
^c{E) = h'<^c{e), *ph(£^) = 6$ph(£). (24) 

Equations (22) can be written in a more transparent 
form if we introduce the effective energy dependent elec- 
tron relaxation time: 



r(e) 



elmujB 



Me), 



where Mc is given by Eq. (6). Then 



'a 

P \ = I h(e-A/)/T 

X) Jrnc-\{e-p)yT 

K(e)T{ey_ (_dl 
de 



(25) 



(26) 



where the energy integration represents the statistical av- 
eraging of the relaxation time with the energy derivative 
of the Fermi-Dirac distribution. Note that Eqs. (22) and 
(26) are valid for any electron degeneracy. Note also that 
Eqs. (19) and (41) of Paper I should contain A/'c(e) instead 

of Tie. 

Generally, r(£) is an oscillating function of electron en- 
ergy as discussed in Sect. 4. Let us mention one important 
case when T{e) and Mc{e) vary with e much slower than 
df{e)/de. Then, for strongly degenerate electrons, one ob- 
tains the results which look formally similar to those for 
B = 0: 



e^c^ricTip) 
a K, , K 

a K, — In 

3e de 



A 



^Vt 



3e2 
AAe(e)T(e) 



(27) 



with rie — J^o{p)- 

The equations of this section are equally valid for the 
quantizing and non-quantizing magnetic fields. 

4. Relaxation time T{e) (or ^{E)) 
4-.1. Quantum oscillations 

Evaluation of the transport coefficients cr, a, and k con- 
sists of two stages. First, the equations of Paper I for 
</3j7nss'(e) should be solved and the function <I>(e) (or, 
equivalently, ^[E), or T(e)) determined. Second, the en- 
ergy integrations (22 ) have to be performed, and the trans- 
port coefficients ( |l9| ) found. The second stage corresponds 
actually to the energy averaging of the relaxation time; it 
will be analysed in Sect. 5. Here we consider the ffi'st stage. 
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The main problem is to reproduce correctly quantum 
oscillations of "^(E) or T(e) which occur since electrons 
of energy e populate new Landau levels with growing e. 
Population of an n-th level takes place when the energy 
variable z/ given by Eq. ^ exceeds n. The oscillations orig- 
inate from the square root singularities of the density of 
states of free electrons in a magnetic field. Just behind the 
threshold for a given level (0 < (z/ — n) ^ 1) <&(£) and T{e) 
behave as ^/v — n which makes the oscillations important 
especially for not too high n. 

Below we obtain practical equations for '^{E) or T(e) 
at various electron energies. 

4-2. Semiclassical approach for v ^ 1 

Itv ^ 1 (electrons occupy many Landau levels), ^{E) can 
be calculated numerically, using the technique of Paper I, 
but it cannot be expressed in a closed analytic form. We 
will derive (Sect. 4.4) accurate fitting expressions based 
on the semiclassical approach described below. 

When many Landau levels are occupied, the relaxation 
time r(e) is expected to be close to the familiar classical 
non-magnetic quantity 



^(e) = UiVaatrie), 



(28) 



where vq is the electron velocity without magnetic quan- 
tization effects {vq = poc^ /e, po — rncPo — mcy/ i?^ — 1), 
rii is the number density of scatterers, and crti (e) is the 
transport scattering cross section: 



o-tr(e) 



J ^ J df]'f7(e)(i-cose). 



(29) 



In this case a{Q) is the differential scattering cross sec- 
tion, 8 is the scattering angle, dfi and dfl' are the solid 
angle elements of non-quantized electron momenta p and 
p' before and after scattering, respectively. In the Born 
approximation. 



a(e) 



\U{q)\ 



2g2 



1 



^0 • 2 



e 



(30) 



Assuming the isotropic distribution over momentum 
transfers hq = p' — p in the non-magnetic case, Eq. ( p9[ ) 
can be rewritten as 



crtr(e) = J 



a{Q) 



An 



(31) 



Now consider weakly quantizing magnetic fields. Using 
Eq. (4) , we can introduce the quantized transverse momen- 
tum (n) — pq sin ?9 — ^2mhijjBn, where n is a Landau 
number, and d is an electron pitch angle. Accordingly, ev- 
ery n corresponds to two values of i?, below and above 7r/2. 
Then an integration over ^ from to tt can be replaced 
as: 



dd sva-d 



mfiujB 
Po\Pz\ 



n{E) 



E E 

n=0 r;=±l 



(32) 



At this stage we replace the integrals over and -d' in 
Eq. ( pl| ) by the sums over n, rj and n', ry'. The dips of 
T{e) or <&(e) behind every Landau threshold are caused 
by the density-of-state singularity when either n or n' 
equals n{E). Since n{E) is assumed to be large, these 
terms correspond to i? « tt/2 or « 7r/2, respectively. 
The sum over rj or rj' in these terms is equivalent to in- 
troducing a factor 2. Thus, the term n = n{E) can be 
evaluated by setting i? = 7r/2 in Eq. (^l|), multiplying by 
2mhuJB / iPo\Pz\) and integrating over all remaining vari- 
ables. The term n' — n{E) is similar. Then 



where 



crci(£) = J dfl' cr {e){l^ cos e), 

CTq(£) = [ di7' CT(e) cos^ I?' = / dn'aie) 



sin^ e 



(33) 

(34) 

(35) 



The last term in Eq. ( |33| ) comes from the terms n = n{E), 
n' = n{E). Its denominator contains the longitudinal mo- 
mentum pz at n — n{E) which vanishes just at a new Lan- 
dau threshold and produces the required quantum oscilla- 
tions. The first term, crci(e), is non-oscillating. Actually it 
should be somewhat lower since we have subtracted the os- 
cillating term from ctr but we neglect this difference in the 
present section. Equation ( ^ ) can also be inaccurate for 
small momenta transfers associated with the transitions 
between discrete neighboring n and n' (which require a 
more detailed consideration). In order to allow for this 
effect, we introduce the lower momentum transfer cutoff 
<Z = <Zmin while calculating i7ci(e) and crq(£). We specify 
(7min in Sect. 4.4. 

From Eq. (O), for the scattering on ions we obtain 



'^d."q(e) = 47r (^-^j R%{E,y), 
where 

i?S(i?,y) = Ai-4a2, 

i?^(iJ,y) = y)_A2 + ^A3, 



(36) 



A. ^ - 



j'' dz 



[z + uf 



Ai = - In 

2 \y + u 

A2 = u m 



2{l +u){y + tiy 



A, - 



1-r 



3 2, 

—u In 

2 \y + u 



y + uj ■ 2{1 + u){y + u)' 

- u{l - y) 
l + u\ u^{l-y) 



2(l + u)(y + w)' 



(37) 
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u = [hqs / {2pq)]'^ is the screening parameter, and y = 

[hqmin/i2po)?. 

For the scattering on impurities, fT™''(e) and (T™P(e) 



are obtained from Eq. by replacing Z — > (Zi,„p — Z) 
In the case of the scattering on phonons, from Eqs. ( p^ 
and (13^) we obtain 



(38) 





8n ( Ze^\'^ o „h 


where 






— r r 


Rt{E,y) 


= <(i?,y)-iL2 + 


Lk 


= fc / z'^-^e-^'^dy, 


Li 


w 




= ^ [(l + u;y)e— ^- 


L3 


3 

= ^ [(2 + 2wy + 
- (2 + 2w + u;2)e-'"] 



^0 

3c2 



(39) 



w = (4/3)(po7'T2/?i)^- As explained in Sect. 3.1, these 
equations are strictly valid T-q < T < Tm, and pro- 
duce reasonable estimates at T^/ < T <C Td . 

Now, using Eqs. (|2^ ) and (|3^), we can easily express 
the inverse energy dependent relaxation time as a sum of 
the main (non-oscillating) and oscillating terms. 

J^.S. Non-magnetic case 

If _B = the relaxation time of an electron is readily 
obtained from Eqs. (^), (|3|) by omitting the oscillating 
term and setting qmin — (jj — 0). Then we immedi- 
ately recover Eqs. (23) - (28) of Paper I for T(e) due to 
the Coulomb and phonon scatterings. 

The function 'i'{E) is defined by Eq. (^). The function 
Mcie) which enters Eq. ( |25| ) and relates r(e) and ^'(i?) is 
given by Eq. (6). If i? ^ 0, A/'o(£) is easily calculated by 
replacing the sum over n with the integral, which yields 
the trivial result Afde) — (po/Ti)^ / (Stt'^). The appropriate 
functions '^{E) are 



p6 
^0 



3A{E)E^ ' 



*ph(i?) 



pi 
^0 



3£;2i?p^(£;,o)' 



(40) 



where A{E) = R^^(E,0) is the Coulomb logarithm given 
by Eq. (24) of Paper! 

Once the energy dependent relaxation time To(e) is 
known, the transport properties are easily evaluated from 
Eq. (27). In this way we reproduce familiar transport co- 
efficients of degenerate electrons for the scattering mech- 
anisms of study. 



4-. 4- Fitting expressions for v > \ 

We have performed extensive calculations of "^{E) for the 
Coulomb and phonon scattering potentials, using the for- 
malism of Paper I. We have mainly used the distribu- 
tion function framework with the fixed-spin basis. Ac- 
cording to Paper I, this yields sufficient accuracy in most 
cases of interest. However, the accuracy has been addi- 
tionally controlled using the density matrix formalism. We 
have considered a wide range of magnetic fields 10^° G ^ 
B <j 10^'^ G sufficient for applications. In the case of the 
Coulomb potential, ^^q^E) depends also on the screen- 
ing wavenumber (Sect. 3.1). Instead of qs, it is conve- 
nient to introduce the dimensionless screening parameter 
u — \hqs/{2pQ)Y, defined in Eq. (37). We have calculated 
'^g{E) treating u as a free parameter varied from to 
1. These results are equally valid for the scattering on 
ions and impurities (Sect. 3). In the case of the scattering 
on phonons, ^ph(£') depends on the Debye-Waller factor 
defined in Eq. (15). The effect of this factor is described 



by the dimensionless parameter w = (4/3)(po''T2/?i)^ in 
accordance with Eq. (39). While calculating \E'ph(-E), we 
treat u; as a free parameter varied in a wide interval, 
< u; ^ 10. 

Our computations cover an extended range of electron 
energies. The energy variable v has been varied up to 100, 
allowing population of up to 100 Landau energy levels. 

We have produced vast tables of "^{E) (for different E^ 
h, u, w). They are inconvenient for practical use but we 
have been able to fit all the results by analytic formulae. 
We have started with simple semiclassical Eqs. (|2^) and 
(^3|). They reproduce the main features of ^(£') at large ly 
but they appear to be inexact near the Landau thresholds. 
We have managed to modify Eqs. ( ^ ) and (|3^) in such a 
way that the resulting formulae fit '^{E) accurately for 
any parameters &, u, w and for v > 1. We have treated 
the minimum momentum cutoff Qmin, or y (Sect. 4.2), as 
a fit parameter. The fit formula reads (i^ > 1) 



*a(-B) 



^o{E) 



Rci{E,y) 



(41) 

(42) 
(43) 



where y = 0.5(P„/Po)2, P„ = y/E'^ - 1 - 2bn; n = 
n{E) = Int(z^) in accordance with Eq. (^; the functions 
Rd{E,y), Rq{E,y) are defined in Eqs. (37) and (39). 

The expression for ^b{E) depends on n. When n > 3, 
for all scattering mechanisms of study we have 



"fbiE) 



PoPn'fo{E)/i3b) 



[R^{E,y+) + R^{E,y-)]{l + i.{P„/Vbf/^y 



(44) 
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where y± = 0.5[(P„-i ± P„)/Po]^- When n = 1 or 2, we 
obtain 



PoPn 



Q2(e++6,o,i) + Q2(e- + Cs,o,i) 
WTWW+WiE + iW' 



(45) 



where kc — '2 and fcph = 1 according to Eq. (gj), ^± 



(Pq ± P„)V(2&), 4^ 



,g,)V2, = 0, and 



\JficJ[eS) is the magnetic quantum length. The functions 
Q2{^,n',n) have been introduced in Paper I. We have 



0^(^,0,1) 



Ce 



-C 



iC + 



/o C + C 



,dC = (l + Oe«i?i(0-l, 
dC 



9 



(46) 



where g = 1 + w/ (4z^), and Pi(^) is the integral exponent 
which is easily calculated (Abramowitz & Stegun 1972). 

Equations (41) - (42) are valid for the Coulomb and 
phonon scatterings. The function ^'a(£') is associated with 
the modified semiclassical expression of ^(E) (Sect. 4.2), 
whereas vE'b(£') is introduced to improve the fit accuracy 
just behind a new Landau level threshold, at (;/ — n) <C 1. 
In the limit of P — > the equations reproduce the correct 
non-magnetic expression ^(P) = ^'o(P)/Pci(P, 0). 

The fitting equations describe the tables of ^'(P) for 
all values of b, u, w (see above) at any energy P. They are 
based on the correct semiclassical expressions (Sect. 4.2), 
and their accuracy does not become worse with increasing 
P as demonstrated in Figs. || and The figures present 
$(P) versus for different scattering mechanisms, mag- 
netic fields b, and charge numbers Z. For illustration, the 
screening wavenumber Qs in the case of the Coulomb scat- 
tering has been calculated assuming strong degeneracy, 
with the Fermi level being equal to the current value of 
electron energy, ji — e. This allows us to find the elec- 
tron number density from Eq. (5), and determine then qs 
from Eqs. (|l^) and (p^. For the scattering on phonons, we 
have considered several values of T which enter the Debye- 
Waller factor defined in Eq. (|l||). If Ty < 1, the Debye- 
Waller factor is negligible, i.e., = 1 in Eq. (|l^). Fig- 

ure ^ demonstrates population of the low-lying Landau 
levels n=l, 2, 3 and 4. Figure^ displays population of 
higher levels, n = 20 and n = 99. The overall accuracy of 
the fits is around 3% which seems to be satisfactory for 
fitting oscillating functions. The accuracy becomes worse 
(up to 8% for n > 5) at < - n < 1, where *(P) con- 
tains dips, but this is not important: the dips are easily 
smeared out due to subsequent energy averaging of ^'(P) 
(Sect. 5). 

Previously the function $(P) was calculated and fit- 
ted, accurately and thoroughly, by Hernquist (1984). The 



100 



10 



0.1 



0.01 
^ 100 



10 



0.1 



0.01 




Fig. 2. Calculated and fitted function '&(P) versus v for 
Coulomb and phonon scatterings at different parameters of 
matter in the case when the electrons populate the Landau 
levels n=l,2,3,4. 



author considered the Coulomb and high-P phonon scat- 
terings in iron matter at P = 10^°, 3x 10^°, 10", 3x 10", 
10^2, 3 X 10^2, 1013, 3 X 1013, 1014 G and included pop- 
ulation of n < 30 Landau levels. His fit formula (his Eq. 
(178)) contained 3 adjustable parameters for every Landau 
level (90 fit parameters for every B and each scattering 
mechanism). However, Hernquist (1994) did not take into 
account the Debye-Waller factor which made his results 
for the phonon scattering not very accurate (Paper I). Our 
fitting formulae require no tables of fit parameters; they 
are valid for any magnetic field and Landau level number 
of practical interest, as well as for any chemical element. 



4-. 5. Low energies, v < 1 

li V < 1, the electrons with energy e occupy the ground 
Landau level n=0. In this case ^'(P) is expressed in a 
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Fig. 3. Same as in Fig. ^ for population of the Landau levels 
n=20 and 99. 



closed form using the results of Paper I. For the Coulomb 
and phonon scatterings, we obtain 



p2 



2Q2(?++es,0,0)^ 



(47) 



where ^+ = 2Pq/6, and is defined in (45). The function 
Q2(C,0,0) is (cf. Eq. (46)) 



?2(C,0,0) 



Qf(e,o,o) 



1 



ic+cr^'^ € 







^dC - 7-e«i?i(0, 
^dC = e«i?i(gO- (48) 



The electron transport properties in the ultra-quantum 
limit <C 1 are discussed, for instance, by Yakovlev (1980, 
1984). 

5. Energy averaging of relEixation time or \P{E) 

The results of Sects. 4.4 and 4.5 fully determine *(i?). The 
next step is to evaluate the longitudinal electron transport 



coefficients from Eqs. ( p_9D and (22). This should be done 
by a numerical energy integration in Eq. (22). The inte- 
gration corresponds to statistical energy averaging of the 
relaxation time which broadens the quantum oscillations. 
The broadening is produced by the thermal effects: all elec- 
trons with energies |£ — A*| ^ k^T are known to contribute 
into the transport properties. 

Actually, quantum oscillations can also be broadened 
by other mechanisms which have been neglected in Sect. 3 
and 4 (since we have used the Born approximation and as- 
sumed strictly elastic electron scattering, see Paper I). Ex- 
act theory of the broadening would have been very sophis- 
ticated, and it is still absent. We incorporate the effects 
of the broadening in an approximate manner. In addition 
to the thermal broadening, we will take into account the 
collisional broadening of the Landau levels, and the broad- 
ening due to weak inelasticity of electron collisions. 

To include these effects we replace (e — ^) by (e — 
/i)T/rfe in the integrands of Eq. (22), with = r + [(7 + 
Ae)/ (277^3)]- Here 7 is the collisional width of the Landau 
energy levels which we set equal to 7 = Ti/ro, tq being the 
field-free electron relaxation time (Sect. 4.3; Paper I). We 
include 7 into Tk in a way familiar to the semi-quantitative 
treatment of magnetic oscillations in terrestrial metals 
(e.g., Shoenberg 1984). In sohd matter, T < Tm, we set 
7 = 7ph + 7imp (even for T < Ty , ignoring the failure of 
our low-T phonon treatment for very low T, Sect. 3.1). 
The quantity Ae is a typical energy transfer of an elec- 
tron in a collision event. We set (A£)ion ~ 'hiu-p 
(A£)ph ~ hwpGoit), and {Ae)imp « 0^ 
ion plasma frequency, vi is the ion thermal velocity. Go (t) 
(with t = T/Tp) describes reduction of typical frequencies 
of phonons excited or absorbed by electrons at T ^ Td- 
We set G'o(i) ^ G{t), where G{t) is the reduction factor 
introduced in Sect. 3.1. Our choice of 7 and Ae is rather 
phenomenological. However the results (the transport co- 
efficients) are not too sensitive to this choice. If more ac- 
curate values of 7 and Ae appear in the future, they could 
easily be incorporated in our calculation scheme. 

Using Eqs. (22), and we obtain the foUow- 



where Up is the 



ing practical equations for the longitudinal electric con- 
ductivity CT, thermal conductivity k and thermopower a 
as functions of electron chemical potential /x, temperature 
T, magnetic field B, and nuclear composition of matter: 



a ^ Iq, k = 



3e2 



'J K,^ ^ 'J ex. 1 

e 



(49) 



J. 



ion 
^0 



imp 



(eC - 1)2 



A 



2.473 X 10^2 s- 



(■^imp -^)^ ^imp 
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1.794 X 10^^ — s-^ 



(50) 



In this case C, — {e — fJ-)/ikBTk), and 9 = k-BTk/{mc^). 

Equations (49) allow us to evaluate the required trans- 
port coefficients for any parameters of study indicated in 
Sect. 2 and 4.4. 

The integration in Ij has been optimized taking into 
account oscillating behaviour of '^{E) and the square root 
singularities at the Landau thresholds En = \/l + 26ri. 
The integrals are presented as sums of separate integrals 
over all those intervals [£'„, £",1+1] which fall in range from 
E = E^ - 4061 to = + 4061, with E,, = n/{mc^). 

For high enough temperatures, when this range cov- 
ers more than 25 thresholds of the Landau levels n, in- 
tegration in Ij over every selected interval [£"„, is 
carried out using Eq. 25.4.34 of Abramowitz & Stegun 
(1972) with 4 mesh points between neighboring Landau 
thresholds. This procedure is specifically adopted for inte- 
grand functions which behave as \/E — En in the vicinity 
oiE = En. 

For low temperatures, a more accurate numerical inte- 
gration is required near E — E^^. For this purpose, we se- 
lect those intervals [£'„, En+i] which span the range from 
E = E^ — W to E^ + W. One or two of the selected in- 
tervals which include the points E — E^ ± 86 are further 
subdivided by these boundary points. We integrate over 
any such interval using the 128-point Simpson formula. 
The integration over intervals not spanning the selected 
range is done as in the case of high temperatures. 

This algorithm ensures fast and accurate (with error 
^ 0.1%) evaluation of the integrals Ij. 

In order to combine the scattering on low-temperature 
phonons and impurities at low T we ignore the breakdown 
of the low-T phonon approximation at T ~ Tu. We eval- 
uate separately the integrals Ij {j—Q, 1, 2) in Eqs. (50) 
due to the phonon and impurity scatterings at T < T^. 
Then we calculate the total values of Ij for T < as 



-1 _ /tp^\ 



.J {Ij )~ + {Ij )~ , and determine the transport 
coefficients from Eqs. (49) and (50). 

For practical purpose, we should also calculate the 
electron number density Ue (which specifies mass density 
p = munc{A/Z)) and its derivative with respect to the 
chemical potential that determines the electron screen- 
ing wavenumber Qc according to Eq. (|lj). The practical 
expressions for Uc and due/dn are given in Appendix C 
of Paper I. In these expressions we will also include the 
broadening of the Landau levels. The thermal broaden- 
ing is naturally implanted in these equations just as in 
Eq. (22). The inelastic scattering can be important for the 
kinetic properties but is inadequate for thermodynamics. 
Thus, we incorporate the coUisional broadening by replac- 
ing T -> Tth = T + [7/(27rA:B)] in Eqs. (CI) and (C3) of 
Paper I. 

In practice, one often needs the transport coefficients 
as a function of mass density p rather than the electron 



chemical potential p,. The inverse dependence of p on p 
and T can be found by iterations from the fitting for- 
mula for ne{fJ-) given in Paper I. For example, if the elec- 
trons populate several Landau levels and p > mc^, we 
use the following rapidly converging algorithm. The zero- 
order value of p can be set equal to pq = C\J {nic)'^ + ppg, 
the field-free value (^) in the limit of strong degeneracy. 
At any subsequent iteration step i = 1, 2, 3, . . ., we calcu- 
late ne'"* from the fitting formula with p — pi-i. Then we 
determine pi for the next iteration as pi — Ca/ {mc)^ + pf, 

taking po = Pfo and pi ^ pi^i{ne/ne~^^)^^^. 
6. Results and discussion 

Using the results of Sect. 4 and 5, we have created a com- 
puter code which calculates the longitudinal electric and 
thermal conductivities and the thermopower of degenerate 
electrons for densities 10^ g cm~^ <p^4xl0^^g cm~^, 
magnetic fields B ^ 10^^ G, and arbitrary nuclear com- 
position in the outer crust of a neutron star (as described 
in Sects. 2, 3.1, and 4.4). The code is based on theoretical 
formalism that is strictly justified for the Coulomb scat- 
tering of electrons on ions at T > T^, for the scattering 
on high-temperature phonons at Td ^ T ^ Tm, and for 
the scattering on impurities at T < Tm. The results are 
also valid (semi-quantitatively) for the low-temperature 
phonon scattering at Tu ^ T <C Td as discussed in 
Sect. 3.1. 

Figures || - § show typical density dependence of the 
longitudinal transport coefficients in a neutron star crust 
with a strong magnetic field. 

Figure^ demonstrates quantum oscillations of the elec- 
tric and thermal conductivities with increasing density 
in '^^Fe matter for B = 10^^ G at two temperatures, 
T ~ 10^ K and 10^ K. Every oscillation is associated with 
population of a new Landau level, starting from the first 
excited level n = 1 a.t p ^ pB ^ 5x 10^ g cm~^ (cf. Fig. 1). 
In the displayed density range (10^ g cm~^ < p < 3 x 10^ 
g cm~^), up to about 60 Landau levels are populated. For 
comparison, dashed lines show the conductivities CTcI and 
Kci in the absence of the magnetic field. 

For the lower temperature, T = 10^ we have T ^ Tb 
(see Fig. 1). The adopted magnetic field B = 10^^ G is 
quantizing (Sect. 2), and the oscillations are quite pro- 
nounced. At this temperature, matter is solid (T < Tm), 
and the oscillations are produced by the scattering of elec- 
trons on phonons. These oscillations are significantly am- 
plified by the Debye-Waller suppression of the electron 
scattering (Paper I) . In the domain of strongly quantizing 
magnetic field {p ^ pB, Sect. 2), the transport coefficients 
differ drastically from those at i? = 0. With increasing 
p in the domain of weakly quantizing fields {p 3> pb-, 
T <C Tb): the oscillations become weaker, and the con- 
ductivities are seen to be close to the non-magnetic ones. 
Every oscillation of the electric conductivity contains one 
maximum and subsequent dip, while the oscillation of the 
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thermal conductivity may show a secondary maximum as 
explained, for instance, by Yakovlev (1980, 1984). The 
amplitudes of quantum oscillations of the electron con- 
ductivity are always stronger than those of the thermal 
conductivity. 

For the higher temperature T — 10^ K in Fig. ^, the 
magnetic field B = lO^'^ G is less quantizing (the ratio 
Tb jT remains larger than unity but it is 10 times smaller 
than for T = 10^ K). Accordingly, the quantum oscilla- 
tions are noticeably weaker due to the thermal broadening 
of the Landau levels. The broadening is more pronounced 
in the thermal conductivity than in the electric conduc- 
tivity. The jumps of the transport coefficients at p « 10^ 
g cm~^ are associated with solidification. For lower p (at 
given T), matter is melted, and the electron relaxation is 
produced by the Coulomb scattering on ions. For higher 
p, ions solidify, and the electrons scatter on phonons just 
as at r = 10^ K. At higher T ^ Tb, the oscillations would 
be entirely smeared out by the thermal broadening. Note 
that the transport coefficients in the ion liquid with the 
magnetic field at T = 10^ K are systematically larger than 
those in the non-magnetic matter. This means that the 
transport coefficients averaged over magnetic oscillations 
are larger than the coefficients for B = 0. This differ- 
ence between the oscillation-averaged and non-magnetic 



coefficients is explained by the classical effect of the elec- 
tron Larmor rotation on the Coulomb logarithm in the 
transport cross section for the electron-ion scattering: the 
non-quantizing magnetic field decreases the Coulomb log- 
arithm (and amplifies the electric and thermal conductivi- 
ties) by reducing the maximum effective impact parameter 
(e.g., Yakovlev 1980). 

Figure || shows quantum oscillations of the transport 
coefficients in iron matter with weaker magnetic field, B = 
10^2 G, at lower T (10*^ and lO'^ K). Since the magnetic 
field is lower, the oscillations associated with population 
of the same Landau levels occur at lower p (cf. Fig. 1). If 
T = 10^ K, matter is solid in the displayed density range, 
10"* g cm-''' < p < 10^ g cm-3, while for T = 10^ K it 
solidifies at p « 8 x 10^ g cm~'^. Qualitative behavior of 
the curves is the same as in Fig. |[ 

Figure ^ presents quantum oscillations of the ther- 
mal and electric conductivities in carbon matter for B = 
10^2 G at T = 10^ and lO'^ K. The main features are the 
same as in Figs. ^ and p|, a lthough the melting tempera- 
ture r,n (r = 172, Eq. (0)) is lower than for iron. Pop- 
ulation of the Landau levels is mainly independent of ion 
species {A, Z), and occurs at the same p as in Fig. ||. If 
T = 10^ K, matter is melted, in the displayed parame- 
ter range. In the case of T = 10^ K, matter is liquid for 
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logp (g cm ) 



Fig. 5. Same as in Fig. for B = 10^^ G at T = 10^ and lO'^ K. 




Fig. 6. Same as in Fig. ^for ^'^C matter. 
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Fig. 7. Dimensionless longitudinal thermopower for ^^Fe and ^'^C matter, B — 10^^ G and T — 10® and 10^ K. Dashes show 
the _B = curves. 



p ^ 10^ g cm^"^ and solid at higher p. In this case quantum 
oscillations are pronounced stronger than at T = 10^ K. 

Figure |^ shows quantum oscillations of the dimension- 
less longitudinal thermopower ae/fce (defined by 
Eq. (49)) for B = 10^^ G. The upper panel corresponds 
to iron matter at T = 10'' K (cf. Fig. Two lower panels 
correspond to carbon matter for the same conditions as in 
Fig. ^. The thermopower is known to be more sensitive 
to the electron scattering mechanism, than the conduc- 
tivities. Its quantum oscillations are seen to be stronger, 
and more complicated. The oscillating thermopower dif- 
fers significantly from that at _B = 0, and it can change 
sign (Yakovlev 1980, 1984). The jumps associated with the 
solidification are pronounced weaker than in the conduc- 
tivities. 

The curves displayed in Figs. |^ - are calculated tak- 
ing into account coUisional and inelastic-scattering broad- 
ening of the Landau levels (Sect. 5). These broadening 
mechanisms appear to be less significant than the ther- 
mal broadenings for the scattering on ions in Coulomb 
liquid and on phonons in the crystalline matter. Never- 
theless non-thermal broadenings become more important 
with decreasing temperature. In particular, the inelastic- 
scattering broadening becomes significant at T <g; Td, 
when the employed approximation of high-temperature 
phonons breaks (Sect. 3.1). The coUisional and inelastic- 



scattering broadenings also become important with de- 
creasing density near the domain of incomplete ionization 
(Sect. 2). 

Figure ^ shows oscillations of the electric and thermal 
conductivities produced by electron scattering on charged 
impurities in iron matter with B = 10^^ G at low tempera- 
ture T = 10^ K assuming, for illustration, that the scatter- 
ing on phonons is inefficient. The impurity number density 
is supposed to be riimp = O.OSni, and the impurity charge 
number is Zimp=28. The impurity screening wavenumber 
gimp in Eq. (|l|) is set equal to gimp = (47rnimp/3)^/^. 
At this low temperature, the coUisional broadening of the 
Landau levels becomes more important. With decreasing 
T, the thermal broadening would die out, and the shapes 
of the oscillations would be entirely determined by the col- 
lisional broadening. Note that we include the non-thermal 
broadening mechanisms in an approximate manner. This 
introduces an uncertainty into the transport coefficients 
at low temperatures. However the case of low tempera- 
tures is not very important for applications. Anyway, the 
number density and charge number of impurities are not 
known, and this introduces much larger uncertainty in our 
knowledge of the transport properties of very cold stellar 
matter. 
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Fig. 8. Same as in Figs. § - ^ for electron scattering on charged impurities in ®^Fe matter at _B = 10^^ G. 



7. Summary 

We have obtained practical formulae for evaluation of the 
longitudinal electric and thermal conductivities and lon- 
gitudinal thermopower in degenerate dense matter (10^ 
g cm~'^ ^ p ^ 4: X 10^^ g cm~^) of outer neutron star 
crusts with strong magnetic fields B = 10^° - 10^'* G 
(at stronger fields, the splitting of electron Landau lev- 
els concerned with anomalous electron magnetic moment 
should be taken into account). The results are expressed 
in terms of the energy averaged electron relaxation time 
r(e) or the function ^'(-E). We have found (Sect. 4) accu- 
rate and simple analytic fits for ^>{E), which are valid in 
wide ranges of the parameters of stellar matter for three 
electron scattering mechanisms: for the Coulomb scatter- 
ing on ions in gaseous or liquid phases, for the scattering 
on high-temperature phonons in solid matter, and for the 
scattering on charged impurities in solids at low temper- 
atures. We have proposed (Sect. 5) an efficient energy av- 
eraging procedure which allows us to evaluate rapidly any 
longitudinal electron transport coefficient. The numerical 
examples are given in Sect. 6. 

Note that our formalism can easily be extended to 
higher densities, p > 4 x 10^^ g cm~'^, in the inner crust 
of a neutron star, where free neutrons appear in matter in 
addition to electrons and atomic nuclei. In this case nu- 
clei can occupy substantial part of volume, and one should 



take into account finite sizes of nuclei by multiplying the 
Fourier image of the scattering potential, |?7(q)p (Sect. 
3.1), by a squared nuclear formfactor (e.g., Itoh et al. 
1984). However, this effect can be shown to be not very 
strong. For instance, in the absence of the magnetic field, 
it reduces the effective electron collision frequencies typi- 
cally by about 20 %, in the inner crusts of neutron stars. 
Accordingly, the present theoretical framework is expected 
to be sufficiently accurate in the inner crusts as well. 

The results of this work are required for studying var- 
ious processes in neutron star crusts (e.g., Yakovlev & 
Kaminker 1994). First of all, we mention thermal evolution 
(cooling) of neutron stars. The outer crust produces ther- 
mal isolation of the stellar interior. The electron thermal 
conductivity is most important for calculating the rela- 
tionship between the surface and interior temperatures of 
the star and for evaluating the distribution of the temper- 
ature over the neutron star surface. The latter distribution 
can be strongly anisotropic owing to anisotropic character 
of the thermal conductivity. The anisotropy of the surface 
temperature leads to modulation of the surface thermal 
radiation due to stellar rotation. The modulation has been 
observed with ROSAT in soft X-ray radiation of several 
neutron stars, in particular, PSR 0656-1-14 (Finley et al. 
1992, Anderson et al. 1993) and Geminga (Halpern & Holt 
1992, Halpern & Ruderman 1993). Correct interpretation 
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of observations requires the knowledge of the distribution 

of the effective temperature over the neutron star surface 
combined with the solution of the radiative transfer prob- 
lem in the neutron star atmosphere (e.g., Pavlov et al. 
1995). Much work has already been done by Hernquist 
(1984), Van Riper (1988, 1991) and Schaaf (1988, 1990) in 
calculating the temperature profiles in neutron star crusts 
with quantizing magnetic fields, and in analysing the cool- 
ing of magnetized neutron stars. Many important and re- 
liable results concerning the temperature profiles have al- 
ready been obtained (particularly, by Van Riper 1988) in 
the approximation of uniform plane-parallel layer with the 
magnetic field normal to the surface. However the overall 
problem of the surface temperature distribution is com- 
plicated and requires further study (as discussed, e.g., by 
Yakovlev & Kaminker 1994). The solution of this prob- 
lem can be based, to some extent, on the above results. 
Secondly, the electric conductivity and thermopower are 
important for understanding the evolution of neutron star 
magnetic fields (e.g., Urpin et al. 1986, 1994). We plan to 
construct the models of outer crusts of neutron stars with 
strong magnetic fields in our subsequent works. 

The Fortran computer code for evaluating the trans- 
port properties based on the above results is distributed 
freely by electronic mail upon request. 
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